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Abstract 

We discuss the dynamics of parton cascades that develop in dense QCD 
matter, and contrast their properties with those of similar cascades of gluon 
radiation in vacuum. We argue that such cascades belong to two distinct 
classes that are characterized respectively by an increasing or a constant (or 
decreasing) branching rate along the cascade. In the former class, of which 
the BDMPS, medium-induced, cascade constitutes a typical example, it takes 
a hnite time to transport a hnite amount of energy to very soft quanta, while 
this time is essentially inhnite in the latter case, to which the DGLAP cascade 
belongs. The medium induced cascade is accompanied by a constant flow of 
energy towards arbitrary soft modes, leading eventually to the accumulation 
of the initial energy of the leading particle at zero energy. It also exhibits 
scaling properties akin to wave turbulence. These properties do not show up 
in the cascade that develops in vacuum. There, the energy accumulates in 
the spectrum at smaller and smaller energy as the cascade develops, but the 
energy never flows all the way down to zero energy. Our analysis suggests 
that the way the energy is shared among the offsprings of a splitting gluon 
has little impact on the qualitative properties of the cascades, provided the 
kernel that governs the splittings is not too singular. 
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1. Introduction 


The strongly collimated cascades of gluons that are radiated by an ener¬ 
getic parton is an essential ingredient of the physics of jets that are produced 
in various high energy collisions, including heavy ion collisions (for a recent 
review see [1] and references therein). This is because the inclusive distri¬ 
bution of the hnal particles is to a large extent determined by that of the 
radiated gluons at the time of hadronization, a feature known as the local 
parton-hadron duality [ 2 ] . In the case of heavy ion collisions, it is expected 
that multiple interactions of the hard partons with the deconhned matter 
present in the hnal state, induce another type of cascade [amiE]. Its elemen¬ 
tary radiation process is described by the BDMPS-Z theory PITlElIHlIiniE]. 
In the ideal case where only medium induced processes are retained, the re¬ 
sulting cascade exhibits particular features that make it very different from 
the cascade of gluons generated by radiation in vacuum, commonly described 
by the DGLAP evolution equation [T2[ |T3l [H]. In particular, the BDMPS 
cascade exhibits scaling properties, akin to wave turbulence na, that are ap¬ 
parently absent in the DGLAP cascade. Also, as was recently emphasized, 
the BDMPS cascade provides a natural mechanism for the transport of en¬ 
ergy towards large angles, which contrasts with the strong angular ordering 
of QGD cascades in vacuum [121 HZl [TSl US] . 

The main motivation of the present paper is to analyze the origin of the 
qualitative differences between the two types of cascades, with the goal of 
fostering our understanding of the BDMPS cascade. This cascade, as we have 
just recalled, may play an important role in determining the properties of jets 
produced in ultra relativistic heavy ion collisions. Insight will be obtained 
by comparing its properties with those of the more familiar DGLAP cascade. 
As we shall see, a quantity that plays a crucial role in this comparison is the 
characteristic time, that it takes a gluon of energy u to split into two 

other gluons. We refer to this time as the branching time. In the BDMPS 
cascade the branching time is an increasing function of cu, G(a;) ~ yG;, so 
that it decreases along the cascade as gluons of smaller and smaller energies 
are emitted. In other words, the branchings accelerate, and as a result it 
takes a hnite time to transport a hnite amount of energy from the leading 
particle to gluons carrying vanishingly small amounts of energy. This time, 
which we shall refer to as the stopping time [ 20 ], is inhnite for the DGLAP 
cascade in which t^{u) is independent of u. 

Further features of the cascades follow from this crucial difference in the 
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branching times. Consider for instance the energy distribntion D{u). In the 
BDMPS cascade this exhibits a persistent scaling behavior all the way down 
to a; = 0, of the form D{u) ~ while in the DGLAP cascade scaling 

cannot be achieved in a hnite time. The emergence of snch a scaling relates 
to the existence of a stationary solntion to the eqnation that governs the 
evolntion of the energy distribntion as a fnnction of time. This stationary 
solntion is of the form 

D,,(u) = Ad, (1,1) 

UJ 

and is associated to a constant (independent of uj) flow of energy, 


^(o;) 


d£:(u;) ^ ujDstj^^) 
dt U (w) 


= const. 


( 1 . 2 ) 


where £ (cn) stands for the total energy carried by glnons with energies larger 
than cn, i.e., £{oj) = f^dLj'D(uj'). The existence of a constant flow of energy 
at arbitrarily small cj implies that 

£(uJo) = [ duDstiuj) = f —t^{uj) (1.3) 

J loq Jojq ^ 

remains hnite when ujq tends to zero. This forces t*(ci;) to be an increasing 
fnnction of uj, which is indeed the case for the BDMPS cascade, for which 
f*(a;) ~ y/u. In the DGLAP cascade t*(a;) is constant, i.e. independent of 
u, and the integral of the stationary scaling spectrnm diverges. As we shall 
see, a scaling spectrnm D{u) ~ l/cn wonld also be expected in this case, bnt 
it takes an inhnite time to develop. 

The properties of the cascades depend also on the splitting kernel, that is, 
on the way the energy is distribnted between the offsprings dnring a splitting. 
However, this turns out to have a minor effect on the main characteristics of 
the cascade, as compared to that of the transport time scale just mentioned, 
at least as long as the splitting kernel is not too singular. In that case 
the cascades develop as if the branching were completely democratic, with 
the two offsprings taking each half the energy of the parent gluon. The 
interactions responsible for the splittings can then be considered as local 
(in energy space). As already mentioned, these properties of the BDMPS 
cascade that we have briefly listed, are typical of wave turbulence [T5] . 

This paper is organized as follows. In the next section, we briefly recall 
the main features of the BDMPS-Z theory, and the properties of the ideal 
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medium induced QCD cascade. By ideal we mean that we ignore the fact that 
the BDMPS mechanism is valid only within a range of energies of the emitting 
gluons, while we shall allow gluons to split all the way down to zero energy 
according to the BDMPS mechanism. This introduces small, unphysical, but 
well identihed effects, that play however no role in the present discussion. We 
contrast the resulting properties of this ideal cascade with that obtained by 
solving (a simplihed version of) the DGLAP evolution equations. In Section 
3, we analyze the two types of cascades from a more general point of view. 
We isolate in the splitting kernel a contribution that we associate with a 
characteristic branching rate, and a function that determines the probability 
that one of the offspring carries a given fraction z of the parent gluon. We 
show that the main distinction between the two cascades comes from the 
branching rate, and in particular whether it increases with decreasing energy 
or not. 

This general analysis is completed in Section 4 by a study of explicit 
models, that can be solved analytically, and that provide insight into the 
emergence of the scaling solutions. In these models we examine the two 
situations of constant or accelerating branchings, hrst in the case of complete 
democratic branching where the offsprings carry a fraction z = 1/2 of the 
parent energy at each splitting, and then in the case where the distribution 
of energy fractions z is independent of z. This analysis conhrms the general 
discussion of Section 4. The main conclusions of the paper are summarized 
in Section 5. The paper contains also four technical appendices. Appendix A 
presents the construction of the analytic solution of the BDMPS cascade both 
in the case of an inhnite medium and that of hnite medium, and Appendix 
B presents the solution of the simplihed DGLAP evolution equation that we 
use in the main text. Appendix G and D provide more details on the analytic 
solutions of the models that are discussed in Section 4. 

2. The ideal in-medium QCD cascade 

We consider an energetic (on-shell) gluon, with initial energy E, trav¬ 
eling through a medium, for instance a quark-gluon plasmej^ This gluon, 
referred to in the following as the leading particle, interacts with the plasma 
constituents, and these interactions induce radiation. The radiated gluons 


^Our main results do not depend crucially on the detailed properties of the medium 
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Figure 1: Illustration of a gluon cascade that is initiated by a gluon with energy E. 
Four generations are displayed. The branching time Ai^, that corresponds to the lifetime 
of generation i, decreases after each branching as in a BDMPS cascade. The inclusive 
distribution D{x,t) measures the probability to find in the cascade, at time t, a gluon 


with energy xE. The rate equation (2.5) describes how this distribution evolves with time 
t. 


themselves radiate, eventually generating a cascade. Our goal is to study the 
average properties of such a cascade, focusing in this paper on the inclusive 
gluon distribution, integrated over transverse momentum, 

dA^ 

(2.4) 

dx 

where x = u/E is the energy fraction of the gluon observed at some time t 
along the cascade, cf. Figfor an illustration. {D(x,t) can be viewed as an 
energy density, with D{x,t)dx being the energy contained in mode^with 
energy fraction between x and x + dx.) As was shown in [2T], D{x, t) obeys 


^We use the term “mode” as a synonymous for radiated gluon throughout this paper 
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the following equatioi]|^ 

= im, (2.5) 


where t* is a characteristic time scale (to be discussed shortly). The right 
hand side of Eq. (2.5), I[D], has the form of a ‘gain term’ minus a ‘loss term’. 
The gain term describes the production of a new gluon with energy fraction x 
via the decay of an ancestor gluon having energy fraction x' = x/z > x. The 
loss term describes the disappearance of a gluon with energy fraction x via the 
decay x —)■ {zx , (1 — z)x}, with 0 < 2 ; < 1. In writing the equation above, we 
have used the symmetry of the kernel under the exchange z ^ 1 — z. This has 
the consequence that the original singularity at z = 0 has been transferred 
to z = 1. As written, the loss term is divergent at z = 1. However, it is 
easily checked that this divergence is cancelled by an analogous one present 
in the gain term (see Ref. ini for details). 

The characteristic time 


= -\,(S) = i./?, (2.6) 

a ayq 

where a = asNc/n, with the number of colors, will be referred to as 
the stopping time |2U]. It is the time at which most of the energy of the 
incoming parton has been radiated into soft gluons (with typically x < 0.1). 
This time differs by one inverse power of the coupling constant from the 
typical time, if takes the gluon of energy E to branch into 

two gluons with momentum fractions of order one, in short the time it takes 
the initial parton to make one democratic branching. Let us recall that in the 


•^As originally derived, this equation is written in terms of the light-cone time x"*" = 
{t + z)l\j2 « with t the physical time. Similarly, its original kernel involves the 

component « PoV^ of the four momentum. Then is given by 
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Thus, modulo a redefinition of q (q —>■ qV^), one can reads Eq. (2.51 in terms of the physical 
time t and energy E. This is the convention that we adopt throughout this paper. 
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BDMPS-Z theory, the branching time = -s/uj/q of a glnon with energy 

u results form the matching between the usual formation time r = uj/k\, 
and the typical tranverse momentum squared acquired during that time r 
via collisions with the medium consituents, i.e., = qr, with q the jet 

quenching parameter. 

The parameter q, called the jet-quenching parameter, is a characteris¬ 
tic transport coefficient that enters the BDMPS-Z mechanism. It controls 
the momentum broadening (the average momentum squared acquired by a 
parton propagating in the medium for a time r is {k±)^ ~ qr (as just used 
above), as well as the energy lost by this parton during its propagation, 
AE ~ gr^. It is valid for frequencies, cjbh ^ ^ ^ where (parametrically) 
Uc ~ with L the length of the medium, is the maximum energy that 
can be taken away by a single gluon (the mechanism that dominates the en¬ 
ergy loss). When oj > ojc, the spectrum of medium-induced radiation drops 
rapidly. The lower limit is that of incoherent emissions, and is reached when 
the branching time is of the order of the mean free path between successive 
collisions. For energies uj < cubh, the suppression due to coherence between 


successive emissions (the so-called Landau-Pomeranchuck-Migdal effect) dis¬ 
appears. In this paper, we shall ignore this lower bound, and assume that 
the medium induced gluon splittings is described by the BDMPS spectrum 
all the way down to zero energy, and refer to the corresponding cascade as 
to the ideal BDMPS cascade. 

In writing Eq. (2.5), we have isolated the dependence of the BDMPS ker¬ 
nel on the energy fraction of the parent gluon, the factors 1/y/x in Eq. (2.5), 
from the 2 ; dependence of the splitting process. Thus, the kernel IC{z) in 
Eq. (2.5) is a “reduced” kernel given by 


1^(1 - 2 ) 13/2 | 2 ( 1 - 2 )] 3 / 2 ' 


(2.7) 


This kernel collects contributions from the z dependence of the actual branch¬ 
ing time (left out in the dehnitiorj^ of f*, Eq. (2.6) above), and from the 
leading order splitting function Pgg{z) = Nc[f (z)]"^ / z{l — z) (we restrict our 
discussion to purely gluonic cascades). The factor ^/x may then be combined 
with in order to define an effective time scale for the branching of a gluon 


^The actual branching time for offsprings carrying fractions z and 1 — z of the initial 
energy is -^/xz!! — z)t*. 
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carrying the fraction x of the initial energy, t^{x) = t^y/x. This time scale 
decreases as x decreases, meaning that the rate of emission of soft gluons 
is higher the softer the gluon: it follows that gluon splittings occur faster 
and faster as one moves down the cascade, and, as a result, it takes a hnite 
time (of order f*) to transport a hnite amount of energy form x = 1 down 
to a; = 0 PO]- It is actually convenient to make this x-dependent time scale 


more evident, and write Eq. (2.5) as 


d 


dz }C{z] 


D{x/z,t) D{x,t) 
U{x/z) t^{x) 


dzzlC{z) =X[D]. (2.8) 


This form of the equation plays an important role in the foregoing analysis. 

In the present study, we shall use a simplihed version of the reduced kernel 
lC{z) given in Eq. (2.7), namely 


lC{z) 


^ 3 / 2(1 _^) 3/2 ’ 


(2.9) 


that is, we shall replace the smooth function f{z) in the numerator by 1. The 
main motivation for using this simplihed kernel is that one can then obtain 
the solution analytically. This solution was already discussed in Ref. [5]. 
Its explicit construction is presented in Appendix A (see Eq. (A. 11)). The 
simplihed kernel (2.9) captures properly the singular behaviors at z —)■ 0 
and z ^ 1, while the smooth function f{z) plays a negligible role, as can 
be verihed through a comparison of the analytic solution with the numerical 


solution obtained with the full kernel (2.7). We shall in fact verify throughout 


this work that the main qualitative features of the cascade are to a large 
extent independent of the specihc form of the (reduced) kernel. 

Before we go into a discussion of the main features of this analytic so¬ 
lution, we wish to contrast this cascade of medium induced radiation with 
the cascade of gluons emitted in vacuum by an energetic (off-shell) parton. 
This is described by the DGLAP evolution PI- A simplihed version of the 
corresponding evolution equation for the inclusive one particle distribution 
reads 




dz ^ dz 

D{x/z,t)—a / -- D{x,t), 

Jo ^ — z 


( 2 . 10 ) 


where here the time variable is related to the virtuality of the emitting 
parton, t = IiiQ^/Qq, and again only gluons are taken into account. This 













equation has the form of Eq. (2.8) if one identifies 


U{x) 


= a, 


and 




z(l-z)' 


( 2 . 11 ) 


( 2 . 12 ) 


Thus the DGLAP equation differs from the BDMPS equation in two major 
aspects. First, the kernel JC{z) is less singular near z = 0 and z = 1. Second, 
the rate of successive branchings is independent of the parent energy, i.e., it 
is constant along the cascade. We shall see that the latter property is what 
makes the major difference between the BDMPS and the DGLAP cascades. 


S.I. The ideal medium-induced QCD cascade 

We now return to the ideal medium-induced QGD cascade, and discuss 


the main features of the solution to Eq. (2.5) obtained with the simplified 


kernel (2.9). This solution reads (see Appendix A, Eq. (A.11)) 


D{x,t) = 


T 


\fx (1 — x)T2 


exp 


-TT- 


X 


A 


(2.13) 


This solution exhibits two remarkable features: a peak near x = 1 asso¬ 
ciated with the leading particle, and a scaling behavior in l/\/x at small x 
where the x dependence factorizes from the time dependence, i.e. 


r 


D{x, t) ^ —= e 

'X 


(2.14) 


An illustration of this solution is given in Fig. left panel. The energy 
of the leading particle, initially concentrated in the peak at x < 1, gradually 
disappears into radiated soft gluons, and after a time t ~ (i.e. r ~ l/0r 

0.5) most of the energy is to be found in the form of radiated soft (x < 0.1) 
gluons. This is also the time at which the peak corresponding to the leading 
particle disappears (see Fig. |^. These are the reasons that motivate calling 
G the stopping time. At the same time the occupation of the small x modes 
increases (linearly) with time, keeping the characteristic form of the scaling 
spectrum. When the peak has disappeared, the cascade continues to lower 
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Figure 2: (Color online.) The function D{x,t) (Eq. ( |2.13| )) at various times. Left panel: 
the filling of the modes, which proceeds till the disappearance of the leading particle 
peak. The values of r are, for the thick (blue) curves, from bottom to top: 0.1, 0.2, 0.3 
(during this stage the leading particle acts as a source for soft gluon radiation), and for 
the thin (black) curves, from top to bottom: 0.5, 0.7, 0.9,1.0,1.1,1.2 (the leading parton 
has exhausted its energy and the peak has disappeared, while energy continues to flow to 
small X, the amount of energy in each mode decreasing exponentially fast). Right panel: 
energy is constantly injected into the system by a source located at a: = 1 (see Eq. (2.21)). 
After a transitory regime, characterized by a uniform increase with time of the scaling 
spectrum, the system reaches a steady state. The values of r are, from bottom to top: 
0.1,0.2,0.3, 0.4, 0.5, 0.6,0.7,0.8,0.9,1.0. 


X, causing a uniform, shape conserving, decrease of the occupations, and a 
flow of energy towards small x that we now analyze. 


Energy conservation is explicitly implemented in Eq. (2.5) at the level of 


individual splittings. To see how this conservation law manifests itself more 
globally, it is useful to consider the flux of energy, J^(xo,r), towards values 
of X smaller than a value xq. This is defined by 


7{xo,t) = 


d£{xQ,T) 

dr 


£{xo,t) = / dxD{x,T), 


(2.15) 


fxo 


where S{xo,t) is the amount of energy contained in the modes with x > 
Xq, and J^{xo,t) is counted positively for energy moving to values of x < 
Xo (hence the minus sign in the definition of J^. These quantities can be 
calculated explicitly. We have for instance 


£{xo,t) 


dxD{x,T) = e erfc 


^ XO 


/ / TTXq 

vVi-^o 



(2.16) 
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(2.17) 


with erfc(a:) the complementary error function: 

erfc(a;) = / due““ . 

ix 

We note that the total energy “stored in the spectrum”, namely 

lim £’(xo, t) = (2.18) 

decreases with time, and accordingly there is a non vanishing flux of energy 
reaching x = 0 

= 2nte-^^\ (2.19) 

It follows that the complete, energy conserving, solution involves a contribu¬ 
tion 

Dc{x) = nc{T)6{x) with nc(r) = 1 — (2.20) 


somewhat analogous to a condensate where the radiated energy accumulates 
(see next section). Note that when r ~ 1 /a/^, corresponding to the disap¬ 
pearance of the leading particles into soft radiation, about 60% of the initial 
energy has flown into the condensate. 

It is interesting to consider also the situation where the leading particle 
is replaced by a source that injects energy at a constant rate ^ at x = 1. In 
this case, we are led to look for the solution of the following equation 


^D{x,t) = ^h(l - x) +I[D], 


( 2 . 21 ) 


where I[D] denotes the r.h.s. of Eq. (2.5) (or Eq. (2.8)). The exact solution 
of Eq. (2.21) with initial condition D{x,t = 0) = 0 reads 


Dst{x,T) = 


27r-\/x(W^a)) 




( 2 . 22 ) 


and is plotted in Fig. right panel. As time goes on, this solution converges 
to the stationary solution (^/27r)/■^x(l — x), keeping the shape of the small 
X spectrum, with just an overall time-dependent scaling. Remarkably, the 
scaling form D{x,t) = f{t)/^/x is reached well before the stationary state is 
achieved. When the steady state is reached, all the energy provided by the 
source flows throughout the entire system towards the condensate at x = 0, 
that plays the role of a sink, while the population of the various modes stays 
unchanged. 
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2.2. The DGLAP cascade 


We now contrast these properties with those of the DGLAP cascade. The 

using a Mellin transform. 


solution of Eq. (2.10) is obtained in Appendix B 


The initial condition reads D(z/, 0) = 1, with D{u, 0) the Mellin transform of 
D{x,0) = 5(1 — x). The solution can be expressed as the following integral 
(see Eq. (B.4)) 


D{x, t) 


fc+ioo 


dl/ 


27ri 


: exp 


-(- 0 ( 1 ^) + 7 ) f + z/ln 


1 

X 


(2.23) 



Figure 3: Model of a DGLAP cascade according to Eq. (2.23) (left panel), and 
with a source added at a: = 1 (right panel). The curves corresponds to r = 
0.2, 0.4,0.6, 0.8,1.0,1.2,1.4,1.6,1.8, 2.0 from bottom to top. In the right panel, the emer¬ 
gence of the scaling solution D{x) ~ l/x\s clearly visible, as well as the persistent deviation 
from it of the true solution at very small x. 


This integral is not easy to calculate in general. However, it allows us to 
verify a few properties that are relevant for our discussion. First, it is easy to 


check that the energy is conserved by the evolution. Indeed, from Eq. (B.4) 
we observe that 


T(0,f)= / da;D(a;, t) = D(l, t) = 1. 


(2.24) 


Second, the following asymptotic behavior is derived in Appendix B 


see 
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Eq. (B.8)): 


D{x, t) 




In 1/x 


(In 1/x 3> t). 


(2.25) 


Here, in contrast to what happens with the BDMPS cascade, the time de¬ 
pendence of the energy distribution does not factorize, and no simple scaling 
behavior emerges at small x. The growth of the spectrum at small x is in 


fact tamed by the exponent in Eq. (2.25), making the spectrum integrable 


when In l/x > t: all the energy remains in the spectrum, and no energy flows 
to X = 0, unlike the BDMPS cascade. 

A numerical calculation of the DGLAP cascade is presented in Fig. In 
particular, in the right hand side panel, the model with a source is given. 
The solution in that case approaches the stationary solution D{x) ~ A/x at 
not too small x, but it takes a long time to populate the very small x modes 
according to this 1/x law. In fact, for any time t, the true solution deviates 


from the stationary solution in the same fashion as in Eq. (2.25). 


The properties of the two cascades that we have discussed in this section 
can be understood in a broader context that we now describe. 


3. Further properties and general considerations 

The two cascades described in the previous section are governed by the 


same general equation, Eq. (2.8) with, in each case, a specihc kernel IC{z) 


and a specihc time scale f*(x). In order to get further insight into the generic 
features of these cascades and understand why and how they differ, it is 
useful to extend the discussion and regard JC{z) and t*(x) as independent 
quantities. The kernel }C{z) controls how, in a given splitting, the energy 
is shared between the two offsprings. The time scale f*(x) controls the rate 
at which successive splittings occur, and this is allowed to depend on x. 
The independent choices made for }C{z) and for f*(x) may appear somewhat 
artihcial since in a given held theory the two quantities are determined by the 
same microscopic physics, as the examples treated explicitly in the previous 
section show. This separation, however, allows us to better illustrate the 
specihc roles of JC{z) and of G(x). 

Equation ( 2.8[ ) has a number of interesting properties that we now exam¬ 
ine. By direct integration of Eq. (2.8), assuming that D{x, t) is regular in the 
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interval [0,1], it is easy to check that the total energy is conserved, namely, 




dxD{x, t) = 0. 


(3.26) 


Since energy is continuously moving towards the lower values of x, the dis¬ 
tribution will generically develop a divergent behavior at small x. The effect 
of this divergence, in particular on the regularity of D{x, f) as x —)■ 0, will be 
discussed shortly. 


One can also verify that Eq. (2.8) admits an approximate hxed point 


solution, which we refer to as a scaling solution, of the form 




U{x) 


X 


(3.27) 


Indeed, in the region where D{x,t) ~ Dsc{x,t), there is complete cancella¬ 
tion between gain and loss terms. Note that this compensation occurs only 


ioT z > X in the integrals of Eq. (2.8), so that the scaling solution cannot be 


an exact hxed point of the equation. We shall see however that the general 
solution is driven to this approximate hxed point, and it approaches it in the 
allowed region of small x. 


Depending on the behavior of the scaling spectrum (3.27) when x —)■ 0, 


the approximate hxed point solution may lead to a divergent expression for 
the total energy T(0, t) contained in the modes between 0 and 1. In order for 
this not to happen, either the hxed point solution is never reached at very 
small X, or t*(x) is an increasing function of x. Both cases have already been 
met. The former occurs in the DGLAP cascade for which t^{x) = const., 
while the latter is that of the BDMPS cascade for which f*(x) ~ ^/x. To 
make things more precise, it is convenient to consider the energy S{xQ,t) of 
the modes with x < Xq, with xq small but hnite, and its time derivative, the 


how J^(xo,t) (see Eq. (2.15)). A simple calculation yields 


rd../cw Td^A^, 
A A t.(x) Jo A tJx) 


(3.28) 


When Xq —)■ 0, the second term becomes negligible (assuming that the x- 
integral is not too singular as Xq —)■ 0). The hrst term has a simple interpre¬ 
tation: it is the integral of the loss term D{x, r)/t*(x) over all values of x for 
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which zx < Xq. If Xq is small enough, and in the region where both xq and 
xq/z are in the scaling region so that we can replace D{x,t)/ t^{x) by 1/x, 
the flow is simply given by 


'XQ 


dz zlC{z) In - = X. 

"Z 


(3.29) 


The limit xq —)■ 0 exists provided zlC{z) is bounded from above by 1/z in 
the limit z —)■ 0, the latter condition on the kerne|^ insures the quasi-locality 
of the branchings. For the BDMPS cascade, v = 27r, while for the DGLAP 
cascade, v = 7r^/6. 

The emergence of the scaling solution is best seen when we replace the 
leading particle by a source that injects energy at a constant rate ^ at x = 1. 


That is, we look for the solution of Eq. (2.21). In this case, not only a scaling 


solution emerges, but it evolves into a true fixed point solution, namely a 
stationary solution. When a source of energy is present, the flux of energy 


contains, in addition to the contribution given by the first Eq. (2.15), the 
contribution of the source. That is 


J'(xo,r) = A 


dS(xo,T) 

dr 


I[D], 


(3.30) 


' XQ 


where I[D] denotes the right hand side of Eq. (2.8). The corresponding flow 
of energy remains then given by the same formula as above, Eq. (3.28). For 
the BDMPS cascade, an explicit calculation yields 


J'(xo,r) = A 


1 — e 


erfc 


TTXq 

1 - Xo 


-r 


(3.31) 


This formula reproduces the result already mentioned, namely the existence 
of a flux of energy down to xq = 0. At large times, all the injected energy 


flows to X = 0, which plays the role of a sink. Eq. (3.31) also reveals how the 


Xq dependence gradually disappears as we approach the stationary solution. 


^This property may be considered as a precise definition of what we call quasi- 
democratic branching. In this situation, the role of the splitting kernel is secondary as 
compared to that of the branching time in determining the qualitative properties of the 
cascade. These properties are then qualitatively similar to that of a cascade dominated 
by completely democratic branching with a kernel peaked at z = 1/2, or in other words 
by local (in energy space) interactions 
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This is illustrated in Fig. In the BDMPS cascade, the scaling regime 
extends all the way down to x = 0. This implies the existence of a hnite flow 
of energy at x = O’*" and the formation of a condensate (or a sink) at x = 0 
where the injected energy is accumulating (see Eq. (2.20)). 

When the stationary solution is reached, the flow is equal to that of the 
energy injected in the system, namely A. Equating this with the general 
expression (3.28), and assuming xq to be small enough so that both xq and 
xq/zq are in the scaling region, we fix the normalization of the stationary 
solution in the scaling region. 



Figure 4: The flow in the presence of a source A = 1. Left: as a function of Xq for various 
times (r = 0.01,0.02,0.05,0.1,0.2,0.5,1, from bottom to top). Right: as a function of 
time for various values of xq (xq = 0.1,0.05,0.01,0.001 from left to right). The lowest 
curve represents the limiting behavior of the flow at xq = 0, obtained from Eq. (3.31). 


For the DGLAP cascade, in the presence of a source, one expects a steady 
spectrum of the form 

Af 

Dst{x) = (3.33) 

vx 

However this scaling spectrum is achieved only when In 1/x <C In 1/xmin ~ t. 
For smaller x, the 1/x divergence is tamed, in such a way that the integral 
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giving £^( 0 ,t) remains finite. As a result there is no flow of energy reaching 
X = 0. These features are visible in the plots of Fig. Thus, in DGLAP 
like cascades, there exists a dynamically generated infrared cutoff, ~ Xmin, 
that prevents the total energy to diverge, as it would if the scaling solution 
were to exist all the way down to x = 0 . At the same time, it is necessary 
for the integral to diverge (in the absence of this cutoff) in order to allow 
for an inhnite amount of energy to be stored in the spectrum. Since the 
hnite integral must account for the total energy that has been injected in the 
system after a time t, At, we have 


At 



rsj 


At^ 1 
-In- 

min 


(3.34) 


so that In l/xmm ~ t, or Xmin e“^h Because all the injected energy can be 
stored in the spectrum within the interval [xmin, 1 ], there is no flow at x = 0 . 
The cutoff Xmin decreases rapidly with time, thus allowing more and more 
energy to be stored in soft modes. 

The discussion above leads us to expect generically two types of cascades. 
To keep the discussion simple, we shall assume the following scaling form for 
the characteristic time: t*(x) = t^x'^. The two classes of cascades correspond 
to specihc values of 7 : 

Type A: 7 < 0 

The integral of the scaling solution diverges when x —)■ 0. At any hnite 
time, the scaling spectrum does not extend to arbitrary low x, but only to 
the value Xmin- In the presence of a steady source, the energy £^(xmin, t) that 
is stored into the system is proportional to the time t, 


^(^min; t) ~ t ~ 



rsj 


ln(l/Xn 


x: 


for 7 = 0 
for 7 < 0 


(3.35) 


We can thus deduce the value of Xmin in the two cases: when 7 = 0, Xmin ~ 
e“^* and when 7 < 0 , Xmin ~ either case the point x = 0 is reached 

only in the limit f —)■ cx). 

Type B: 7 > 0 

In this case the scaling spectrum is integrable in the infrared. This rehects 
the fact that it takes only a hnite time to transport a hnite amount of energy 
from X = 1 to X = 0, so that the energy does not accumulate in the spectrum 
but at X = 0. The cascade also involves a characteristic time, or equivalently 
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a characteristic value of x, where multiple branchings become dominant. 
This is given by t ~ ~ t^x'^, so that 


x, = — 


1/7 


(3.36) 


When X 3> only few branchings occur and the spectrum is well approxi¬ 
mated by the lowest order in a perturbative expansion, 


D{x,t) ~ Di{x,t) ~ —x/C(x). 

f * 


(3.37) 


It is in the multiple-branching regime, i.e., for x < x*, that the scaling 
spectrum develops. 


L>sc(a:,t) = ^^^/(i), 


X 


(3.38) 


where f{t) is a function of time that depends on the kernel ]C{z). Note that 
Eq. (3.38) indicates that x^ —)■ 1 as 7 —>■ 00 , reflecting the fact that the 


scaling regime develops on the entire available energy range. 

At this point, it is useful to comment briefly on a peculiarity of the 
BDMPS cascade, that makes the transition between the dilute, single-branching 
regime, and the the multiple branching regime invisible in the spectrum 
(2.13). Indeed, for x^ -C 1 (or equivalently r <C 1), the spectrum is of 


the form 


D{x,t) ~ 

lx 


(3.39) 


and this is so independently of whether x S> x^ or x x. This property 
is specihc to the BDMPS spectrum for which x/C(x) ~ l/\/x, so that the 
leading perturbative result has already the scaling form. 


4. Simple models 

In order to illustrate the above general considerations, we consider in this 
section simple models that differ by the form of the branching time f*(x) 
and by that of the kernel lC{z). We assume, as we did at the end of the last 
section, that t*(x) scales as t*(x) = t*x'>', with 7 an arbitrary parameter in 
the physically interesting range 0 < 7 < 1. We shall consider two models 
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that illustrate the generic types of cascades that we have identihed in the 
previous discussion. In the hrst model, to be referred to as model ^4, 7 = 0, 
corresponding to a constant branching rate, as in the DGLAP cascade. In the 
other model, to be referred to as model B, 'j = 1/2, the branching rate grows 
as in the BDMPS cascade. This difference in the branching rate, coded here 
in the value of 7 , has a profound impact on the solution: if the branching 
time is independent of x, it takes an inhnite time to populate the smallest x 
modes. In the case where 7 > 0, the splitting rate increases as we move down 
in X, and it takes only a hnite time to transfer energy from the highest to the 
lowest values of x. To emphasize that the x-dependence of the branching rate 
is the dominant factor that determines the main features of the cascade, we 
shall consider two extreme versions of splitting kernel. The hrst corresponds 
to democratic branching, with /C(z) peaked at the value 1/2. The second is a 
hat kernel, independent of 2 ;. We start by analyzing models with democratic 
branching. 


4-1. Democratic branching 

Complete democratic branching corresponds to a kernel of the form 


lC{z) = 25{z-l/2), 


dz z'IC{z) = 1 . 


(4.40) 


With this choice of /C, Eq. (2.8) becomes 

ot t^[2x) U[x) 


(4.41) 


where the theta function in front of the gain term reminds us that the support 
of D{x, t) is limited to a: < 1. This theta function will often be omitted in the 
following (the constraint on the support of D being implicitly implemented 
in the function D itself). 

The two models that we consider correspond to diherent choices for t^{x). 
For model A, t^{x) = A is constant, and Eq. (|4.41|) reads 


drD{x, t) = 2D{2x, t) — D{x, r), (model A) (4.42) 


with T = t/t^. It admits the following, non integrable, hxed point solution, 
D{x) = For model B, f*(x) = t^y/x, and Eq. (4.41) reads 


dT-D{x,T) =—= V2D{2 x,t) — D{x,t) 

Dx L 


(model B) (4.43) 
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with again t = t/t^. It admits the fixed point solution D{x) = ^. In 
contrast to model A, the singularity at small x is now integrable. 

When the energy is injected in the system at exactly x = 1 (either by a 
leading particle or by a source), the distribution function has support on a 
discrete set of points. It is of the form 


D{x,t) = ^ Cfc(r)(5 


fc =0 


—;- X 

2k 


(4.44) 


The equation (|4.42|) then leads, for both models, to the following recursion 

(4.45) 


relations for the coefficients Ckir) 

C^ = A- Co, Ck = Ck-ia’^-^ - Cka\ 

witlfl 


model A: a = 1, model B: a = 2'*', (4.46) 


and A the strength of the external source. When A is non zero, the initial 
condition is Ca, = 0 for all k. When ^ = 0, the initial condition is Co = 1 
and Ck = 0 for all k ^ 0. In either case, the equation for Co(t) is the same 
for both models, and it can be solved easily: 


^ = 0: Co(r)=e-", Aj^O: Co(r) = ^(1 - e""). (4.47) 


When ^ = 0, the exponential decay of Co(r) represents the decay of the 
leading particle. When A ^ 0, this decay competes with the feeding of 
the zeroth mode by the source. The fixed point solution corresponds to 
Ck-i = aCk- Whether it is reached or not depends on the initial condition, 
and on 7 . The flow at Xq = 1/2^° is given by 


fco 


^(fco, T)=A-J2CkiT) = a’^°Ck,{T). 


(4.48) 


This equation has a simple interpretation: the flow of energy that moves to 
modes with values of k higher than k^ is equal to the loss term, a^°CkQ, of 
the equation (4.45) for the coefficient Cfep(r). 


We now analyze successively the solutions of models A and B. 


®In this section, we focus on the case 7 = 1/2. 


presented in Appendix C 


A solution for general values of a is 
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4-1.1. Model A 

We first examine the solution corresponding to a leading particle carrying 
initially all the energy. The initial condition in that case is V(x,t = 0) = 
6(1 — x). The coefficients Ck obey a Poisson law, 

C'fc(r) = (4.49) 


The energy of the leading particle (k = 0) decreases exponentially with time, 
Ck=o = as already noticed. The energy of the mode ko goes through a 
maximum; it increases till r = ko as it is hlled by energy coming from the 
modes with k < ko-, and then it decreases exponentially when the energy of 
these feeding modes is exhausted. It takes therefore an inhnite amount of 
time to populate the modes with ko ^ oo. The total energy is conserved, 
Jq dxD(x,T) = 1 = J2T=o^k(T)- the very small x (large k) modes are in 
fact never populated. Near a: = 0 there is no flow, since for any given r, 
Ck{T) —)■ 0 as A; —)■ oo. Note also that it takes a hnite time, r ~ /c to excite 
modes with a given k\ this is reflected in the fact that in Eq. (4.49) giving the 
coefficient Ckir), the hrst k coefficients of the Taylor expansion in powers of 
time vanish. All these features can be recognized in Fig. upper left panel. 


Consider then the case with a source A ^ 0. The C^’s are then given by 


Ck(T) = A 


k 


1 - 



n=0 



(4.50) 


They become eventually all equal to ^ as r —)■ cxd. However reaching this 
hxed point takes an inhnite amount of time. Indeed, the function of r that 
multiplies ^ is a growing function of r that saturates at 1 for large r, the 
time at which saturation is reached being a growing function of k\ the small 
k modes are rapidly hlled, but it takes a time of order k to hll the modes 
with large k. The how at Xq = 1/2^° is equal to a^°Cfcp(r) = Ckoir) and, for 
any hnite r, it vanishes for ko 3> t, because the ko modes have not have time 
to get excited. These features are illustrated in Fig. upper right panel. 


4.1.2. Model B 

The solutions for model B, with or without a source, are obtained in 
. They are plotted in Fig.[^ lower panels. The striking diherence 
with model A resides in the fact that these solutions quickly approach a 


Appendix C 
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Figure 5: The coefRcients Ck for the solution with a leading particle (left panels) and with 
a source (right panels), for the following values of k, 1, 2, 3,4, 5,6, 7, 8,9, from left to right. 
The top panels correspond to model A, the lowest panels to model B. 

scaling regime, where the time dependence factorizes and the shape of the 
spectrum gets frozen: it becomes independent of k for k large enough. This 
is intimately connected to the fact that, in model B, it takes a hnite amount 
of time to transport energy from a; = 1 to a; = 0, in contrast to model A 
where this time is inhnite. 

This difference between the two situations is best illustrated by consid¬ 
ering an even simpler version of the previous models, in which the splittings 
occur at discrete time steps. After n time steps, the energy of the leading 
particles is shared equally by 2"' particles. It takes an inhnite number of time 
steps to transfer all the energy of the leading particle into an inhnite number 
of particles carrying a vanishing amount of energy. However, whether this 
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infinite number of time steps takes an infinite amount of time depends on 
the model considered. In model A, we assume that the time interval be¬ 
tween two splittings is a constant At. In this case, it takes an inhnite time 
T = reach the end of the cascade. In model B, we assume that 

the time between two successive branchings decreases as one moves down the 
cascade, according to the relation 

At = Atq , a = 2^, (4. 



where Atq is the time it takes to make the hrst branching. The time T/v it 
takes to go through N successive branchings is given by m 


Tn = 



Ato 


l-g-^ 
1 — a~^ 


(4.52) 


For model A, = NAtq diverges as N ^ oo. However, for 7 = 1/2, 
limAr^.oo7jv is hnite, and T/v converges in fact rapidly towards its limit. For 
instance, for a = \/2, l — a~^ > 0.9 already for N = 7. This number A^ = 7 is 
roughly the number of curves that one can distinguish in Fig. lower panels. 
Thus in this model B, the entire energy of the leading particle is transferred, 
in a hnite time, to an inhnite collection of particles carrying each a vanishing 
amount of energy. This is the analog, in this discrete model, of the conden¬ 
sate that we discussed earlier: the energy injected in the system eventually 
ends up at a: = 0 in a hnite time. 


The fact that energy is transferred to x = 0 in a hnite time also helps us 
to understand the factorization of the time dependence and the emergence 
of the scaling solution. The plots in Fig. [^indicate clearly that the hlling of 
the modes takes place in an essentially hnite time, after which the solution 
acquires a (nearly) factorized form, with the time dependence of a^Cfc(r) 
becoming independent of k, that is Cfc ~ 1/a^. In the case where energy is 
injected into the system by a source ^ 7 ^ 0 , a stationary regime is eventually 
reached, but the scaling solution emerges well before that, as the plot in 
Fig. [^suggests. 

The scaling behavior observed in Fig. implies that the dominant con¬ 
tribution of the coefficients Ck given by Eq. (C.7) tends to a factorized form 


77/M- 


(4.53) 
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Figure 6: (color online.) The emergence of the scaling behavior for the distribution 
D{x,t) ~ Ck/x as a function of a:: = 1/2*^ for a = (the discrete set of points on 
which the function is defined have been joined by a continuous curve). Left, the lead¬ 
ing particle case: thick (blue) curves are for r = 0.5,1,1.5, 2 from bottom to top. Thin 
(black) curves are for r = 2.5,3, 3.5,4.5, 5, 5.5 from top to bottom. One sees the disap¬ 
pearance of the leading particle and the uniform growth with time of the small x (large 
k) modes. When the leading particle that was feeding energy into the small x modes 
has disappeared, the energy continues to flow to a: = 0, the populations of the various 
small X modes decreasing uniformly with time. Right panel: Energy is injected in the 
system at a constant rate by a source located at x = 1. The various curves correspond to 
T = 0.5,1,1.5, 2,2.5, 3, 3.5,4.5, 5, 5.5, from bottom to top. 


This can indeed be inferred from Eq. (C.7). At late times, the fnnction /(r) 
is given by the term proportional to e“'^: 



m=l m=l 


(4.54) 


the other terms in the series of exponentials becoming rapidly negli¬ 

gible. One can in fact relate this scaling behavior to that expected for the 
corresponding continnnm solntion of Eq. (4.43). To do so, let us focus on the 
small X region where the density is large enough to allow us to replace the 
sum over k by an integral 



dk = 


1 


dx 

X 


(4.55) 
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with X = 1/2^. We have then 


D{x,t) = ^C'fc(r)5(a; - 2 


k=0 


dx' 


T-tt / -^Ck{x')ir)5{x - x') 


In 2 
1 


x‘ 


xln2 


Ck{x) ■ 


Now, using the fact that a = 2'^, we rewrite Eq. (4.53) as 


so that Eq. (4.56) takes the form 

L)(x,t) 


1-7 


X 


fM 


(4.56) 


(4.57) 


(4.58) 


This is indeed the expected scaling solution (cf. Eq. (3.27)). The continuous 
distribution D(x, r) is plotted in Fig. where the scaling behavior is clearly- 
seen. As we have shown in the previous section, this scaling regime sets in 
for X ^ Xs{t) ~ r^, with Xs separating the single and multiple branching 
regimes (for r ^ 1). To see how this scale emerges in the present model, 
let us recall that the coefficients Ck take the form of a series in with 
a = y/2 > 1, that is 


Ck = '^ B{m, k) < 


(4.59) 


m=0 


where the coefficients B{m,k) can be easily inferred from Eq. (C.7). In 
the continuum limit the sum over the modes m turns into an integral over 
x' = 1/2™, 


1 dx^ T 

~ / —B[m{x'),k{x)]e~^. 


(4.60) 


For a given r ^ 1, the dilue regime corresponds to the region where ^/x < r: 
there one can expand the exponential. Such an expansion is no longer possible 
when when x < Xs{t): then the argument of the exponential is of order unity 
for all values of x' such that a/o/ < r. 
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4-2. Constant kernel 

The previous models capture the main qualitative features of the cascades 
that we have introduced in the previous sections. In order to see the influence 
of the kernel }C{z) on these results, we now examine the case of a flat kernel 

IC{z) = 1, (4.61) 

for which analytic solutions can also be obtained. We shall still consider the 
two cases of model A in which f*(x) = is constant, and model B where 
t*(x) = t^y/x. As before, we set r = t/t^. As we shall see, the qualitative 
features of the cascades of both models are not much altered by the fact 
that the flat kernel allows at each splitting energy to be transported all the 
way down to zero. This is in line with the remark made in Sect. which 
suggests that only a sufficiently singular kernel would lead to a substantial 
modihcation of the cascade properties. 


4-2.1. Model A 


The solution is obtained in Appendix D.l and reads 


D(x, t) = X 



(4.62) 


where Ii(x) is a modihed Bessel function of the hrst kind. The exponential 
decay is slower than in the strict democratic branching case: this reflects 
the fact that very asymmetric branchings are now possible, and these leave 
a large fraction of the energy in the vicinity of the leading particle, which as 
a result is seen to decay at a smaller rate. At large x and not too long time, 
D(x,t) ~ XT, a linear behavior that is clearly visible in Fig. left panel. 
One also observes on this hgure the slow build up of the scaling solution 
(visible around x ~ 0.1). Another important feature of the solution is the 
generation of a dynamical infrared cutoff. This is visible on Fig. (right 
panel), and can be inferred from the asymptotic behavior of the distribution 


(cf. Eq. (B.8)), 


when a; -C e Eq. (D.7), quite analogous to that of the DGLAP cascade 


D(x,t) 


X 


--1 


(4.63) 


This behavior of C(x, r) as x —)■ 0 is what keeps the integral of the spectrum 
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Figure 7: Model A with a flat kernel. Left panel: the function D{x,t) in Eq. (4.62), with 
T increasing from bottom to top, r = 0.1,1,2, 3,4, 5,6, 7,8, 9,10. Right panel: the same 
function D{x,t) in a linear plot, exhibiting the dynamically generated infrared cutoff, in 
agreement with Eq. (4.63). 


finite. As time goes on, the energy injected in the system remains in the 
spectrum, the flow vanishing at a: = 0 


see 


Appendix D.l). 


The same property holds when a source Al5(l — x) is added. In this case, 
the solution converges slowly towards a stationary solution 


A A 

Dstix) = —, 
X 


(4.64) 


where we used the fact that n = 1/4 for a flat kernel (see Eq. (3.29)). As 


already emphasized, this solution is reached only for values of x above the 
dynamical infrared cutoff x^mn (see Eq. (3.35) and the discussion at the end 
of Sect. 1^. 

Ji^.S. Model B 


The solution corresponding to model B is obtained in Appendix D.2 It 
can be written in the form 


D{x,t) = g{x,T) + e ^'^^(l —x). 


(4.65) 


where 

9 {x,t) = /-i/ 2 (r) ^ + /o(r) + /i/ 2 (r) y/x + /i(r)a;, (4.66) 

V X 
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Figure 8: Model B with a flat Kernel. The left panel shows the evolution with time of the 
solution (4.651, with r = 0.1,0.5,1,2,3,4,5,6,7,8,9,10 from bottom to top. The right 
panel shows how the total energy (4.69) contained in the spectrum decreases in time. 


and the various /^(t) are given in Appendix D.2 


behaves very much like that of model A, 


At early times, this solution 


D{x, t) ~ XT. 


(4.67) 


However, as can be seen on Fig. right panel, a scaling behavior in 1/y/x 
quickly develops at small x. The explicit calculation of [Appendix D.2 indi¬ 
cates that the scaling solution is of the form 


D{x,t) ~ 


192 ^ 


(4.68) 


The typical value of x at which the transition occurs is therefore Xg ~ (see 
the discussion at the end of Sect. |^. 

The integral of the spectrum reads: 


dxD{x,T) = 2/_i/2(r) + /o(r) + ^/i/ 2 (r) + (4-69) 

In contrast to model A, this is a decreasing function of time, as can be seen 
in Fig. right panel. Since the energy is conserved at each step of the 
cascade, this implies that the energy flows all the way down to x = 0 where a 












condensate develops. The solution in the presence of a source is determined 
Appendix D.2 The stationary solution in this case reads 


in 


Dst — 


4A 


X 


(4.70) 


and apart from a different x dependence, is similar to that of model A. 


5. Conclusions 

The present analysis clarifies the distinct properties that characterize the 
medium-induced, or BDMPS cascade from the vacuum DGLAP cascade. In 
the first case, the leading parton is essentially on-shell, and radiation is in¬ 
duced by the multiple collisions that the leading parton and its offsprings 
undergo with the medium constituents. In the second case, the leading par- 
ton is off-shell, and the evolution parameter is the virtuality of the emitted 
radiation. We have seen that the two types of cascades may be viewed, for¬ 
mally, in a broader context where the key quantity that differentiates them 
appears to be a suitably dehned branching time, i.e., the typical time it 
takes a gluon to branch quasi-democratically into two gluons with compa¬ 
rable fractions of the initial energy. In the DGLAP cascade, this time is 
constant, while it decreases along the BDMPS cascade. As a result, it takes 
a hnite time to transport a hnite amount of energy from the leading particle 
to very soft gluons, while the corresponding time is infinite in the DGLAP 
cascade. This property is at the origin of the scaling behavior of the energy 
distribution in the BDMPS cascade, which exhibits indeed many properties 
that are reminiscent of wave turbulence. The DGLAP cascade can in fact 
be viewed as a limiting case of the BDMPS cascade, but since the relevant 
transport time scale is inhnite, it takes an inhnite amount of time to reach 
the scaling solution, which is therefore not observed. 

In this paper, we have focused on the energy flow, but the present study 
has a broader scope. In particular, as shown recently, the same characteristic 
features of the BDMPS cascade as exhibited here play an essential role in de¬ 
termining the angular distribution of the produced radiation [mdHKIS]. We 
have discussed only the ideal case where the BDMPS mechanism is allowed 
to work all the way down to zero energy. This is of course an idealization. 
In an actual cascade propagating in a quark-gluon plasma, the energy would 
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not accumulate at zero energy, but rather at the typical scale where dissi¬ 
pation would occur, namely the temperature of the system, that also plays 
the role of the Bethe-Heitler momentum cutoff. While the generic features 
of the BDMPS cascade that we have discussed here are, to a large extent, 
independent of this lower cutoff, the question of how precisely the energy of 
such a cascade dissipates its energy in the medium is an interesting one that 
deserves further studies (for a recent discussion see [23]). 
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Appendix A. Exact solutions: BDMPS-like cascade 


In order to solve Eq. (2.5), it is convenient to perform a change of variable 
that gives the gain and loss terms a more symmetrical form. In the gain term, 
we set ^ = x/z and in the loss term, ^ = xz. A simple calculation shows 
that, in terms of these new variables, the equation reads 


drD{x,T)= d^P{x,OD{^,T) - D{x,t) d^P{^,x), 


(A.l) 


with 


P(x, 0 ^ . (A.2) 

Note that the singularities at z = 1 in the original formulation have been 
moved to .^ = x in both the gain and the loss term. These singularities cancel 
each other when both terms are added. 

In the following we shall consider Erst the case of a large medium, L > t^, 
i.e., Uc < E, and next the case Uc <C E. 


Appendix A.l. The exact solution for Xc > 1 

We consider first the case of a large medium, L > t^., i.e., ujc > E oi 
Xc = oJc/E > 1. In order to calculate separately the gain and the loss terms, 
we may introduce a regularized splitting function Pc{x,f) 

Pe{x,0 = \ - j- -37^- (A.3) 

V X — X + eyi^ 
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Let us analyze the loss term with this regularized splitting function. The 
integral over ^ can be calculated explicitly 



1 

e + e)3/2 


A/e X + e s/e x 


(A.4) 


Note that apart from the fist term that is divergent when e —)■ 0 all subleading 
terms vanish (for any hnite x). Hence, in these variables, the sole effect of 
the loss term is to remove the singularity of the gain term. Accordingly we 
shall write the loss term symbolically as 

r r°° dz 

I (A.5) 

In order to proceed we now introduce the following rescaled distribution 
F(|/, r) = y/x D{x, r), with y = 1 — x. This has the advantage of removing 
the expected singular behavior of D{x) at small x. We get then 


A 

Wr 


F{y,T) 



1 

{y - ifF 


-F{y,T) 


r de 1 

Jy 

(A.6) 


At this point, it is convenient to extend the domain where F{y, r) is dehned, 
from [0,1] to [0, cxd], which amounts to extend the dehnition of D{x) to 
unphysical (negative) values. The motivation for doing so is that it allows 
us to take a Laplace transform. Dehning the Laplace transform as 

poo 

F{a,T)= dye~‘^y F{y,T), (A.7) 

Jo 

we get, for the gain term. 




d( 

(y - 


Fitr) 



dfe-»«FK,T) 



(A.8) 


In this manipulation, the singularity at ^ = y has been moved to z = 0, 
where it takes exactly the same form as that of the loss term, Eq. (A.5), 
with which it therefore cancels when gain and loss terms are added. We then 
arrive at the hnite equation for the Laplace transform F{a,T) 


^F(a,r) = 


dz 


(e — l) F(q!, r) = —2 v^to F(q!, r). (A.9) 
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For a single gluon, the initial conditi on is F{y,T = 0) = 6{y), which yields 
F{a,T = 0) = 1. After integrating (A.9) over t and inverting the Laplace 
transfom, the solution of (A.6) reads 

pc+ioo 1 / 2 \ 


It follows that, in the physical domain, x = l — j/, 0<x<l, 

-2 


D{x,t) = 


T 


^/x (1 — x)^F 


exp —IT 


T 


1 — X 


(A.ll) 


Appendix A.2. Solution for Xc -C 1 

When Xc ^ I, the emission processes are sensitive to the upper cutoff 
on the BDMPSZ spectrum (recall that this spectrum falls very rapidly with 
the energy u of the radiated gluon when u > Uc = qt^)- We implement this 
suppression with a sharp cutoff at a; = Uc, and modify the splitting function 
accordingly as follows 


F(a;,0-^P(x,0^c(x,e), 


6c{x, f) = 6{xc - x) + 6{xc - ^ + x), 

(A.12) 


and the same for P{f,x). It is also convenient to divide the distribution 
D{x,t) into a soft component, Ds{x,t) = 9{xc — x)D{x,t), and a hard 
component, D\^{x,t) = 6{x — Xc)D{x,t). As we shall see shortly, when 
Xc -C 1, the hard component obeys a closed, homogeneous, equation with 
initial condition Dh{x,T = 0) = 5(1 — x), while the soft component obeys 
an inhomogeneous equation where the hard component plays the role of a 
source. Note that the two components are separated by a gap, x G [xc, 1—Xc], 
where D{x, r) = 0. 

Let us consider hrst the hard component, that is the region 1—Xc < x < 1. 
In this region, there is no constraint on the gain term, since f — x < Xc- There 
is a constraint on the loss term however. We write this term as follows: 

An elementary calculation yields 


rx—Xc 


P{^,x) = 


X — 2Xr 


'T* .. / .. / 'T* 'Y* 

JU \/ %jU \/ JU JU f 


lx 


1 Xr 




(A. 13) 


where in the last step, we have used Xc -C 1. In this limit, the correction to 
the loss term is independent of x and the solution Zlh(x,t) reads simply 


Fh(x,r) 


r 


(1 - x)3/2 


TFT 

1 —rc 


(A. 14) 
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Turning now to the soft region {x < Xc), we note that the cutoff affects 
there only the gain term. We write 

dM^,T) = [' d^P{x,ODy,{^,T)+ r\^P{x,OD,{^,T) 

J 1 — Xc ^ X 

diP{i,x)D,{x,T). (A.15) 

The hrst term on the right hand side plays the role of a source term. It can 
be calculated explicitly: 

S{x,t) = -^erfc g2T/\A^_ (A. 16) 



By rescaling the variables, setting ^ = ux^ x = zxc, r = ^/x^Tc, one can put 
the equation for Ds in the form 


dr,Ds{z,Tc) = ^erfc (^/nTc) 
\Jz 


e^""" + 


duP{z,u)Dj^u,T^ 
duP{u,z)P)J^z,T^. (A.17) 


This is the original equation (see Eq. (A.l)) with a source term. Again, 
this can be solved via Laplace transform, using manipulations similar to 
those described earlier. We only write here the hnal solution in terms of the 
original variables. 


Ds{x,t) 



dr'erfc 


e^’^ erfc 




(A.18) 


Note that in the limit a; -C Xc and for r = Tmax = one can expand 

around r = 0. We get 


T>s(x,r) 



1 


a - o? + 0{a^) 

3 


(A.19) 


One recovers corrections that can be obtained by a perturbative calculation 
[22] . More generally, the structure of the cascade is that of a peak localized 
near x = 1, feeding soft modes at small x < Xc- The main features of the 
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resulting cascade of soft modes are quite similar to those of the general cas¬ 
cade that we have detailed earlier in this section. 


We note hnally, that although we have focused on the cascade generated 
by a gluon, a similar analysis can be made for the case where the leading 
particle is a quark. The corresponding energy distribution is obtained by 
using the following simplihed form for the kernel 

^ ^ 1 / 2(1 _ ^) 3/2 ’ (^- 20 ) 

It is easy to check that the quark distribution is related to the gluon distri¬ 
bution as follows Dg{x,t) = xDg{x,t), since the loss term stays unchanged. 
Hence, 


Appendix B. Solution of the DGLAP cascade 

In this Appendix, we write the solution of the DGLAP evolution equation, 
Eq. (2.10) in Mellin space. We dehne the Mellin transform of the energy 
distribution, and its inverse, respectively as 


D{iy,t) = / dxx^ ^D{x,t), 


dz/ 

D{x,t)= —x’'D{iy,t), (B.l) 


where the integration contour in the inverse transform runs parallel to the 
imaginary axis, to the right of the poles of the integrand. By taking the 
Mellin transform of Eq. (2.10), and recognizing the integral representation of 
the Digamma function. 


Uf 


-d. - 7 = 


(B.2) 


one obtains 


^D(z/, t) = - {ij{u) + 7) D{u, t), 


(B.3) 
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where 7 is the Euler constant. It is straightforward to solve Eq. (B.3), with 
the initial condition 0) = 1, the Mellin transform of D{x, 0) =5(1 — 0 ;). 
The solution reads 


D(i/, () = exp [- (V'(i') + 7)*l ■ 


(B.4) 


Transforming back to x space, we write the solution of the DGLAP equation 
as 

r‘C+200 

+ 7 )t + z/In — . (B.5) 

X 


, r™ dz/ 


In general it is difficult to perform the integration over z/ in Eq. (2.23). But 
in some limits one can perform analytic evaluations which are be sufficient 
for the purpose of this work. This is the case in particular for the small x 
limit where one can use the saddle-point approximation. The argument of 
the exponential has a minimum when 

In 1 /x 


^'(g) = 


t 


When In I 3> t, i>(y) ~ —v ^ and '0'(z^) ~ z^ so that z/^ = 


(B. 6 ) 


Hence 


we get 

D{x,t) ~ exp 


— 7 f + 2 y t In — 

X 


rc-\-ioo 


du 


27ri 


: exp 




t 


1/4 


exp 


In^ 1 /x 

Focusing on the leading behavior, we get 


— 7 t + 2\lt In — 

X 


(B.7) 


D{x, t) 




lx 


(B.8) 


This is the standard Double Logarithmic limit of DGLAP. Note that the t 
and X dependence of the energy distribution do not factorize. Obviously, 
since In 1/x S> t, the spectrum is integrable at small x. This can be verihed 
by an explicit calculation. From Eq. (B.4) we have indeed 


E= dxD{x,t) = D{l,t) = 1. 


(B.9) 
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Appendix C. Exact solutions: discrete models with democratic 
branching 


In this Appendix, we provide details on the analytical solution of the 
the discrete models that are introduced in the main text (Section 4). Recall 
that in these models, the distribution D{x,t) has support on a discrete set of 
points, X = 1/2^, the energy in each mode k being denoted by Cfc(r). We set 
here = 1 so that t = t. The coefficients Ck satisfy the recursion relations 


(4.45), that we reproduce here for convenience 

Ck = Ck-ia^-^ 


Co = A-a 


0 , 


Cka^ {k > 1 ), a = 2 ^. 


(C.l) 


In the main text, we introduced two models, referred to respectively as model 
A and model B. In model A, a = 1 (7 = 0 ), which corresponds to a constant 
splitting rate along the cascade. In model B, a > 1 (7 > 0), which implies 
an acceleration of the splittings along the cascade. 

The initial condition depends on whether we consider the system with 
a source, A 7 ^ 0 , or without a source, A = 0. In the first case, the initial 
condition is (^^(t = 0 ) = 0 for all k, while in the latter case, Ck>i(t = 0 ) = 0 , 
C'o(t = 0) = 1, with Co representing the energy of the leading particle. We 
also assume that at large time Ckit) vanishes or is bounded (we verify a 


posteriori that this is indeed the case). We then solve the equations (C.l) 
with a Laplace transform. We dehne, with a slight abuse of notation. 


Ck{uj)= dte-^^Ck{t), Ck{t) = 


'■+*00 1 

— e-^Ckioj). (C.2) 
27r 


Let us hrst analyze the case of a leading particle carrying initially all the 
energy, i.e., the case A = 0. By taking the Laplace transform of Eq. (C.l) 
we obtain 


(l + n;)C'o(a;) = l, uCk{uj) = a’^-^Ck-i{uj) - a^Ck{u;) (fc > 1 ), 

(C.3) 

where we have used Cfc(oo) = 0 and C'fc(O) = dko- The recursion relation can 
be solved easily, with the result 

^ -| 

Ck(cj) = TT — - Co(u), Co(u) = - -. (C.4) 

^ ^ J-i ^ + am uv UV y ^ V ) 

m=l 
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The inverse Laplace transform then yields 


Ck{t) 


r-ioo gCjt 


l-ioo 1 + 1 ^ 



a 


m —1 


Lo + a^‘ 


(C.5) 


For model A, where a = 1, the integrand in Eq. (C.5) has a multiple pole 
of order k + 1 at a; = —1. By closing the contour on the left half-plane {t > 0) 
and picking up the contribution of the residue, one obtains easily 


Ck{t) = ^f"e-k (C.6) 

kl 

In the case of model B, where a > 1, the integrand exhibits a series of 
distinct poles, u = —1, —a, —a^, ■ ■ ■ , —a*'. By closing the contour in the left 
half plane, and collecting the contributions of the residues one then gets. 


Co{t) 

Ck>i{t) 



7n=l m=l 


.m— 1 


-d^t 


1 - 




Note that the exponential decay of Co(f), which represents the decay of the 
leading particle, is the same in both models. The way the population of the 
other modes [k > 0) varies with time depends however on the model in a 
signihcant way. 


The Poisson distribution (C.6) may be viewed as the limit of the general 
distribution (C.7) when a —)■ 1. In this limit, all poles at —a, —a^, • • • ,a^ 
merge together with the pole at a; = — 1, producing a pole of order fc -|- 1 at 
uj = — 1 . 

The other interesting limit is that of large a, corresponding to a rapidly 
accelerating cascade. In this case, the numerical calculation suggests that the 
coefficients Ckit) with k > 1 quickly reach the scaling behavior Ck ~ 1/a^. 
This is illustrated in Fig. |C.9 right panel. In fact when a 3> 1, and for 
not too small times, a^Ckif) ~ aC'i(f), and aC'i(t) itself becomes closer and 
closer to Co(f) as a grows, as can be inferred from the explicit expressions 


Ci(t) = 




^—at 


Co{t) = e-h 


(C.8) 
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At very small t {at <C 1), however, the behavior is reminiscent of that found 
for a = 1, 


Ck{t) = a 


^ Mk-l)/2 i_ _ k{k-l)/2 


1 — 




k\ 


1 — a (/c + 1)! 


(C.9) 


The vanishing of the hrst k coefficients of the Taylor expansion of the coef¬ 
ficient Cfc(r) may be seen as a non perturbative effect, corresponding to the 
fact that modes with high values of k are not excited instantaneously. The 
delay stays however hnite a k grows, in contrast to what happens in model A, 
where it grows linearly with k. These behaviors of the coefficients Ck in the 
two models can be see in Fig. [^of the main text. A similar analysis can be 
done for the case with an external source of strength A. Only the equation 
for the coefficient Cq depends explicitly on A: 


Co{t)=A{l-e-^), 


CAu) =A{- ^ 

00 


1+00 


(C.IO) 


As compared to the case without source, Eq. (C.4), CA^) contains an addi¬ 
tional pole at 00 = 0 corresponding to the stationary source. Using the same 
method as that used in the absence of source, one gets, for model B, 


Ck>i{t) - —j: 

A o'* 


n 

m=l 


■^m —1 


E 




1 ^ a{a^ — 1 ) 

m=l ^ ' 


n 

n^m 


-.n— 1 


(C.ll) 


For model A (a = 1), one recovers by this method the expression (4.50) of 
the main text. As was the case in the absence of the source, when a S> 1, 
for A; > 1 behaves essentially as aCi{t), and eventually as Co{t) for very 
large a and not too small t, with 


CAt)=A 


^—at 




^—at 


(C.12) 


and C'o(t) given by Eq. (C.IO). At small t one observes the same non pertur¬ 
bative behavior as in the case without the source, with Ck{t) given by 


1 ^k{k-l)/2^k+l ^k+2 ^ _ ^k+1 

^ {k + l)\ {k + l)\ 1-a 




(C.13) 


Again, this behavior is visible only in a very small time interval, ta <+i 1. 
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Figure C.9: The coefficients a^Ck{t) for increasing values of a, with a = 1.1,2,5,20, and 
various values of fc (A: = 0, • • • , 10 from left to right). These values of k correspond to 
distinct curves only in the case a = 1.1, where the behavior is close to that of model A 
(a = 1), the energy in each mode reaching a maximum at time t k. The flat behavior at 
small time that develops as k grows clearly shows the delay in the excitation of modes with 
high k values. As a grows, the curves corresponding to the various values of k gradually 
merge together, the merging occurring from the largest k values. Thus for instance for 
a = 2, the curves corresponding to fc > 4 are indistinguishable from that representing 
0^(74(t), for 0 = 5, only the curves representing aCi{t) and 0^(72(t) are visible. As a 
increases further, all curves merge onto that representing aCi{t), the latter becoming 
itself closer and closer to Co(t) = e“*, except at very small t. 
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Appendix D. Exact solutions: constant kernel 

In this section we provide the solntions of the two models A and B with 
a flat kernel 


lC{z) = 1. 


(D.l) 


Appendix D.l. Model A 

The evolntion equation for model A (7 = 0, a = 1) reads 


d 

Wt 




^^D{x, t)=x I t) - ^D{x, t). 


(D.2) 


In order to solve this equation, we perform hrst the following change of 
variables, 

2 /= In - , j/'= In ^ and D(x,t) = xe~^^g(y,t)- (D.3) 

X 4 


Equation (D.2) then yields 


d 


-^9{yA) = j dyg{y,t), 

which can be solved by Laplace transform. We obtain 


(D.4) 


aiy.t) = r = ./I/i(2V(5), 

J_ioo 27 r* V y 


(D.5) 


where Ii{x) is the modified Bessel function of the hrst kind. In terms of the 
original variables, the energy distribution in the cascade reads 


t 


D{x,t) = h I 2y fin - ) e 2 . 


X 


(D.6) 


When X e *, the behavior of D{x,t) is analogous to that of the DGLAP 
cascade (cf. Eq. (B.8)), 


D{x,t) ^ 



(D.7) 
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In particular, it vanishes at a; = 0, that is, an infrared cutoff is dynamically 
generated. This is enough to ensure that the energy of the cascade remains 
in the spectrum, as the following calculation shows 


dxD{x,t) = 


JO 

— t iy- 

2Tii ' ^ 


e-5* / dxx I 

2m 


{u-2)yyr(\/v-\/2)t 


>0 


du 1 
27ri V — 2 


^{l/v-l/2)t 


^ ^(2/u-l)t 


= 1 . 


v=2 


(D.8) 


In the problem with a source A5{1 — x), the function g{y,t) becomes 

which yields the steady state solution 

AA 

Dst = —. (D.IO) 

X 

Note that this solution exists only for not too small x. At very small x, the 
dynamically generated cutoff suppresses the flow of energy, and D{x,t) goes 
to zero. 


Appendix D.2. Model B 

For model B (7 = 1/2, a = v^), the transport equation reads 

^D(a;, 

This is to be solved with the initial condition Dq{x) =5(1 — 0 ;). To solve this 
equation, we start by shifting the distribution by the solution at the singular 
point a; = 1, 

D{x) = g{x) + — x). (D.12) 
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Eq. (D.ll) then yields 


9 , , _i. 

■^ 9 {x,t) =xe 2 +x 


1 


{ 5 / 2 - 


2 ^ 5 - 


(D.13) 


If we use the ansatz g[x) ~ in the r.h.s. of Eq. (D.13), we obtain 


X 


-P- 

2^/x 


X = 


2a — 2> 


2x- 


(D.14) 


Imagine now solving Eq. (D.13) iteratively. The leading order term is ~ x, 
so that in the hrst iteration a = 1, which yields a term proportional to x and 
one proportional to y/x. The next iteration will generate the powers x, \fx 
and P, the third iteration, x, \/x, l/\/x and P. In the fourth iteration, four 
powers are generated, a = —1/2, 0,1/2,1, but the a = —1/2 contribution 
cancels with the second term in the r.h.s of Eq. (D.14) that generates a 1/x 
term. It follows that the iterative procedure does not generate other powers 
than a = —1/2, 0,1/2,1. As a consequence, the solution takes the general 
form 


t) = f-ll2{t) + /o(t) + fl/2{t) Vx + fi{t) X. 
\ X 


(D.16) 


Plugging Eq. (D.15) into Eq. (D.13) and matching the coefficients of the four 
different powers of x we obtain a closed system of equations. 


A 

— e 2 * — 2/i — /i /2 — -/o — -^f-i/ 2 , 


A/2 

= pu 


A 

= 


/-1/2 

0 

(D.16) 

; is convenient to perform a Laplace transform 



POO 

fn{l2) = / dte"^Vn(t). 

Jo 

(D,17) 
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Using the boundary conditions 

fn{oo) = 0 and /„( 0 ) = 0 , 

we then easily obtain 

^/l = 1 ~ 2/i — fi/2 — -fo — 7^/-l/2, 

ly + ^ S 2 ' 


(D.18) 


^/l/2 — 

^/o = 

^f-1/2 = g/o- 


(D.19) 


It is straightforward to solve this system of equations. We hrst get for /i (for 

vyO) 




(D.20) 


which yields 


+ 2z/^ + -z/^ + -u + — ] fi{u) — 


I' + l 


(D.21) 


It is easy to check that the polynomial in the l.h.s of the latter equation has 
a single zero, p = — 1 / 2 , so that 

z/^ + 2z/^ + -z/2 + -z/ + —= fz/+-) . (D.22) 

2 216V 2y ^ ^ 

We can then readily invert the Laplace transform, and get 


flit) = 


dz/ z/^ 


o rco 

Mt) = 2 / 


'-*oo 2 vri (p + lY 
dz/ z /2 


e‘'^ = ij { 24f - 18t^ + 1 e- 2 ‘ 

4! 


27ri 


n5 


-too "' (^ + I) 


. . N 3 P dz/ z/ 
/oW = 7 / 


- .oo 27ri 

. 2 ^ 1 r dp 


n5 


(^+ i) 


n5 


(^+ i) 


e"* = 


e"* = 


e"* = 


2x4! 

3 

4x4! 


i2r - 4r + -r e-^h 
4 


4^3 _ 


1 J.4 -it 

re 2*. 


8x4! 


(D.23) 
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For the case with a source, the equation reads 


^D{x,t) = AS{1 - x) + X [ 


-.D{x,t)- (D.24) 


Following the same steps as before, we look for a solution in the form 
D{x, t) = g{x, t) + (^1 — e“^*j 5(1 — x). 

The coefficient f- 1/2 is given by 


(D.25) 


^ ^ > 2Tii 




4 

A 
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n5 


a{^ + 1) {^ + \) 


16 - ( 16 + 8 f + 2t^ + I e 


The remaining coefficients are then easily obtained 

/i(f) = A(-t^ - -\ -T 

^ ^ V2 6 96 


e 2 


Mt) = 


4 \ „-kt 


e 2 . 


A 


Mt) = ^te 




(D.26) 


(D.27) 


At late times, the constant contribution in the term /_i/ 2 (t) dominates, and 
the solution becomes the stationary solution 




4A 


(D.28) 
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